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ABSTRACT 


We report the detection and mass measurement of a binary lens OGLE-2015- 
BLG-1285La,b, with the more massive component having Mi > 1.35 Mq (80% 
probability). A main-sequence star in this mass range is ruled out by limits 
on blue light, meaning that a primary in this mass range must be a neutron 
star or black hole. The system has a projected separation r_i_ = 6.1 it 0.4 AU 
and lies in the Galactic bulge. These measurements are based on the “microlens 
parallax” effect, i.e., comparing the microlensing light curve as seen from Spitzer, 
which lay at 1.25 AU projected from Earth, to the light curves from four ground- 
based surveys, three in the optical and one in the near infrared. Future adaptive 
optics imaging of the companion by 30m class telescopes will yield a much more 
accurate measurement of the primary mass. This discovery both opens the path 
and dehnes the challenges to detecting and characterizing black holes and neutron 
stars in wide binaries, with either dark or luminous companions. In particular, we 
discuss lessons that can be applied to future Spitzer and Kepler K2 microlensing 
parallax observations. 
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Subject headings: gravitational lensing: micro - binaries: general - stars: neutron 
Galaxy: bulge - black hole physics 


1. Introduction 


All known stellar-mass black holes (BHs) are in close binary systems. Presumably this 
is a selection effect induced by the fact that these BHs have been detected via X-ray emission 
generated by accretion from a close companion. In the near future, the Gaia satellite may 
detect (or place interesting limits upon) BHs in wider binary systems with main-sequence 
companions at semi-major axes 0.2 < a/AU < 5. However, gravitational microlensing 
appears to be the only way to systematically study the populations of isolated BHs and BHs 
in well-separated binaries with dark (BH or neutron star (NS)) companions. This is because, 
in the absence of a luminous companion, BHs do not generate photometric signatures (except 
possibly when accreting from the interstellar medium), and gravitatio nal lensing is unique 
in its ability to detect objects based solely on their gravitational held (Einstein 19361) . 


Gould! (j2000al ) estimated that almost 1% of microlensing events observed toward the 


Galactic bulge are due to BH lenses, with another 3% due to NS lenses. However, even 
if this estimate is correct, not a single such lens has yet been dehnitively identihed. The 
reason is quite straightforward: the principal observable in microlensing events, the Einstein 
timescale is a combination of three physical properties of the lens-source system, 

Oe ,2 _ ^ ^ 4G 
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Here 9e is the angular Einstein radius, vTrei = AU(Zljj — Dg ) is the lens-source relative 
parallax, and Pgeo is the lens-source relative proper motion in the Earth frame. While the 
source distance Ds is usually well-known, none of the other variables in these equations are 
routinely measured. Hence, simple parameter counting implies that two other parameters 
must be measured to determine M and Dl- 

Two quantities that can in principle be measured are the angular Einstein radius 6e and 
the microlens parallax tte, whose magnitude is tte = TTr^i / 9e and whose direction is that of 


the lens-source relative motion fsee lGould fc Hornell2013l Figure 1 for a didactic discussion) 


If these two quantit ies can be measured, then the physical parameters can be disentangled 
flGon]dlh992l . 
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For the great majority of microlensing events, neither tte nor is measured. However, 
for the special case of BH lenses, the Einstein timescales ^e = ^E/^^ cc tend to be 

l ong, and th i s gre a tly enhances the prospects to measure tte via Earth’s orbital motion 
fjSmith et ahl 120031: iGouldl 120041) . Hence, the main method that has been recognized for 
identifying BHs and measuring their masses using microlensing has been to attempt to 
measure 6 *e for these long-timescale BH candidates. 

However, this same large Einstein radius 9^ oc diminishes the already tiny chance 
{p = p = 9^,/6 -e) that the lens will transit a sourc e of radius 6*^,, which is the principal method 
by w hich 9^ can be measured for isolated lenses f)Gouldlll994ai IYoo et al.ll2004 iGould fc Yee 

2ni8h . 


More than a decade ago, when microlensing event detections were almost two orders 
of magnitude less frequent than today, three BH candid a tes w ere identihed, MAGHO-98- 
6, MAGHO-96-5, and MAGHO-99-22. Poindexter et ah 1 2005 ) subsequently showed that 
these had, respectively, low, medium, and high probabilities to be BH lenses. However, none 
were either confirmed or rejected as BHs because there was no measurement of 9^- The 
main idea to measure 9-^ for isolated BHs is astrometric microlensing, which takes advantage 
of the fact that the displac ement of the image centroid f rom the true source posi tions is 
directly proportional to 9e f)Walkeiill995l: iHog et al.lll995l: iMivamoto fc Yoshiil 119951) . Such 
measurements are being pursued by several groups although they are challengin g with today’s 
facilities. However, they may become routine in the future with WFIRST flGould fc Yee 


20141) . 


Black holes in relatively wide (~ 1-10 AU) binaries open a second, less explored path 
toward BH microlensing mass measurements. In contrast to isolated lenses, binaries often 
have large caustic structures, which greatly increases the probability that the source will 
transit these structures. Because the magnification formally diverges to infinity at a caustic, 
a caustic crossing permits a measurement of p = 6'*/6 'e, the ratio of the angular source size 
to the angular Einstein radius, since the observed magnihcation will be affected by the hnite 
size of the source. Because 9^ is almost always easily measured flYoo et al.ll2004|) . this yields 
6 *e. Of course, such measurements require that the events actually be monitored during 
the typically brief (few hour) caustic crossings, which in the past has required either good 
luck or very aggressive followup observations. This situation may change with the ramp up 
of modern surveys like OGLE-IV and KMTNet, which monitor wide fields at greater than 
1/hour cadence for 13 and 16 square degrees, respectively. 


To date, however, no such microlensing BH binaries have been discovered. Part of 
the problem has certainly been missed caustic crossings, but another part is that the en¬ 
tire “paradigm” outlined above actually applies mainly to BHs in the Galactic disk. From 
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Equation ([2]), the microlens parallax tte = . Hence, tte tends to be small for BHs. 

If, in addition, the lens is in the Galactic bulge, e.g., TTrei = 0.01 mas, then an M = 5Mq 
BH would have tte = 0.016. Such small parallaxes are difficult to detect, and especially to 
reliably measure, from the ground. Moreover, at a proper motion /i = 4masyr“^ typical of 
bulge lenses, the timescale would be only ^e = 60 days, which is not exceptionally long. Thus 
Galactic bulge BHs, including BH binaries are much more difficult to detect and measure 
from the ground than Galactic disk BHs. 


Space-based microlensing parallax is well placed to meet these challenges. Rather 
than relying on the fairly slowly accelerated motion of a single observer on Earth 


Gould 


19921) . space-based microl ensing directl y compa.res co ntemporaneous lightcurves from two 


well-separated observers (IRefsdall Il966l : iGouldl Il994bl ) . For an Earth-satellite separation 
(projected on the sky) D^, the microlensing parallax is approximately given by 

fo,© to,sat _ 


ATI 

tte = -p—(A t, A/3); 


At = 


A/3 = - ±M| 


0,sat ? 


( 3 ) 


where the subscripts indicate parameters as measured from Earth and the satellite. Here, 
{to,Uo,tE) are the standard point-lens microlensing parameters: time of maximum0, im¬ 
pact parameter, and timescale. For binary microlensing, these form a subset of a larger 
parameterization, but the formula remains valid. For point-lens events, precisions at the 
level afTTp,) < 0.01 required fo r bulge-BH mass measurements have been achieved in prac¬ 
tice (iGalchi Novati et al.ll2ni5h . For binaries, the challenge is greater bec ause of an intrinsic 
asymmetry in the measurement of (Ar, A/3) (iGraff fc Gould et al.l 120021 ): one linear combi¬ 
nation 


Atcc = (At -|- A/3 cot (/))H 


( 4 ) 


can be measured with exquisite precision from the difference in caustic-crossing times At^c 
seen from Earth and the satellite. Here, 0 is the angle between the source trajectory and 
the tangent to the caustic. However, the orthogonal combination is much more difficult to 
measure. We will discuss this challenge in some detail in Sections 13.41 141 [5l and[6l 


Nevertheless, the main difficulty is the availability of such parallax satellites, which 
must be in solar orbit (or orbiting a solar system body that is not itself orbiting Earth), 
capable of reasonably good photometry in the crowded bulge helds, and, of course, allocated 
to microlensing observations. 


Dong et ahl (120071) made the hrst such microlens parallax measurement using the IRAG 
camera on Spitzer for the microlensing event OGLE-2005-SMG-001. This was, in fact, a 


^For binary lenses, this is replaced by time of closest approach to some fiducial point in the lens geometry, 
which is usually not in fact the maximum. 



























7 


binary star lens and moreover the favored interpretation was a BH binary. Unfortunately, 
however, there was no caustic crossing, so this remains a candidate rather than a conhrmed 
detection. 


In 2014, the Director allocated 100 hours of Spitzer time for Galactic bulge observations 


with the specihc aim of determining Spitzer’s viabilitv as a microlensing parallax satellite. 

Based on this successful test ( 

Udalski et ah 

2ni5a: 

Yee et ah 

2ni5a 

Galchi Novati et ah 


2015 ). which included one mass measurement of a binary f Zhu et ah 20151) . 832 hours were 


awarded for 2015, which is a majority of the 38 days that bulge targets are visible from 
Spitzer due to Sun-angle restrictions. 

While the main focus of this program was to determine the Galactic distribution of 
planets, and the main protoco ls for both Syitz er and supporting ground-based observations 
were determined on this basis flYee et al.ll2015bl) . there was also a signihcant effort to monitor 
binaries, exactly because of the possibility of mass measurements. 


Here we report on the mass and distance measurements of OGLE-2015-BLG-1285La,b. 
The mass of the primary indicates that it is most likely a NS or BH. The system lies < 2° 
in projection from the Galactic center and is almost certainly a member of the Galactic 
bulge population. The held was specihcally targeted by OGLE for factor ~ 5 increased 
cadence to enable early alerts that would permit timely Spitzer observations and to increase 
the probability of resolving unexpected caustic crossings. It was further targeted as part 
of a ~ 4deg^ survey by UKIRT and Wise observatories in Hawaii and Israel, respectively, 
in order to both increase phase coverage and, in the former case, take advantage of the 
capacity of IR observations to penetrate the relatively high extinction in these helds. The 
mass measurement of OGLE-2015-BLG-1285 is a specihc product of these specially targeted 
observations. 


In Section [2] we discuss the observations, with emphasis on Spitzer and the above men¬ 
tioned special targeting. In Section [3l we present a microlens model and demonstrate how 
the physical conclusions follow from the light curve features and the source position on the 
color-magnitude diagram (GMD). In Section 01 we show that future proper-motion measure¬ 
ment of the luminous component (s) of the binary lens will yield a decisive mass measurement. 
In Section |5l we discuss the requirements for breaking similar degeneracies in future events 
using microlensing data alone. Finally, in Section [6l we discuss some other future prospects. 



















2. Observations 


2.1. Special Northern Bnlge Fields 


The high extinction toward the Galactic plane in optical surveys roughly splits the 
Galactic bulge into distinct northern and southern regions. Microlensing surveys tradition¬ 
ally concentrate on southern bulge helds close to the Galactic plane because the event rate 
there is high and the extinction is relatively low. This maximizes the number of detected 
events with the high-quality light curves that are required for planet detection. It is of course 
understood that the northern bulge helds, being roughly symmetric with the south, have just 
as many microlensing events. However, prior to the 2015 Spitzer campaign, only one north¬ 
ern bulge held was targeted for high-cadence observations: OGLE-IV BLG611, centered at 
(£, 6) = (0.33,2.82). Although this held has a long heritage going back almost two decades 
to OGLE-11, no systematic study had ever been made as to which northern bulge helds were 
the most prohtable to target. 


In the course of analyzing the 2014 Spitzer campaign, we realized that Spitzer target 
selection was being heavily inhuenced by optically-based microlensing alerts, whose distri¬ 
bution on the sky was strongly impacted both directly and indirectly by the pattern of dust 
extinction. That is, hrst, to the extent that high-extinction helds are surveyed, it is more 
difficult to hnd and monitor events because they are systematically fainter in the optical. 
Second, because of this very fact, these helds tend to be monitored at lower cadence or not 
at all. By contrast, for any event that can he detected in the optical (i.e., Aj < 4) Spitzer is 
essentially unahected by the dust. In 2015, therefore, special ehorts were made to counter 
this bias, which included taking specihc account of the extinction at each event location 
flYee et al.ll2015bl) . In addition to these general measures, we also identihed several northern 
bulge helds for special observations, including an iL-band survey using UKIRT (Hawaii) and 
an /-band survey from Wise (Israel). One of the four OGLE helds containing these regions 
was BLG611, which was already being observed wit h hourly cadence. Bu t the other three 
helds BLG653, BLG654, and BLG675 (cL, Fig. 15 in Udalski et ah 2015b ') were raised from 
cadences of roughly 0.5/day to 2-3/day. 


The UKIRT/Wise helds were selected using the procedures developed bv IPoleskil (120151) . 
who showed that the product of the surface densities of / < 20 stars and clump stars is a 
good predictor of the microlensing event rate. 
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2.2. OGLE Alert and Observations 


On 2015 June 7 UT 19:39, the Optical Gravitational Lens Experiment (OGLE) alerted 
the community to a new microlensing event OGLE-2015-BLG-1285 based on observations 
with t he 1.4 deg^ cainera on its 1.3m Warsaw Telescope at the Las Gampanas Observatory in 


Ghile ( 

Udalski et ah 

201 

5b|l using its Earlv Warning Svstem fEWSl real-time event detection 

software ( 

Udalski et ah 

1994; 

Udalski 

2003 

). Most observations were in I band, but with 


some V band observations that are, in general, taken for source characterization. These are 
not used in the modeling. At equatorial coordinates (17:39:23.75, —27:49:13.0) and Galactic 
coordinates (0.23, —1.75), this event lies in OGLE held BLG675, with a nominal observing 
cadence of roughly 2-3 times per night. 


2.3. Spitzer Observations 


OGLE-2015-BLG-1285 originally appeared to be a point-le ns event. The proto cols and 
strategies for observing such events with Spitz er are reviewed in IStreet et al.l (l2015l) and are 
discussed in greater detail bv lYee et al.l f)2015bl) . In brief, targets were submitted on Monday 
for observations on the following Thursday through Wednesday for each of the six weeks of 
the Spitzer campaign. 


OGLE-2015-BLG-1285 was selected for Spitzer observations on Monday June 22 UT 
13:33, i.e., for the fourth week of observations, and at standard hourly cadence. Given the 
time of the OGLE alert, it could not have been selected for the hrst two weeks. The Spitzer 
team specihcally considered this event during preparations for the third week but found its 
predicted behavior too ambiguous to select it. Even in the fourth week, it was considered 
highly risky but was chosen specihcally because it lay in a held covered by UKIRT and Wise, 
and so would automatically receive good light curve coverage. Also, it was noted that the 
source was probably a red giant and so about 100 times brighter at 3.6 pm than typical 
targets, even though it was relatively faint in the optical due to high extinction. 


In the sixth (hnal) week, the cadence was increased to 4/day on the grounds that it 
was apparently anomalous. According to the protocols of lYee et al.l fl2015bl) . such increased 
cadence can be used to characterize the anomaly or increase the precision of the parallax 
measurement, provided that the anomaly and parallax are detectable without them. Alto¬ 
gether, Spitzer observed OGLE-2015-BLG-1285 38 times, with each epoch composed of six 
30-second dithered exposures. 
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2.4. Other Survey Observations 


The sky position of OGLE-2015-BLG-1285 was covered by three surveys in addition to 
OGLE, namely UKIRT, Wise, and KMTNet. As with OGLE, the UKIRT and KMTNet 
observations were carried out without consideration of any known microlensing events in the 
held. The Wise observation procedures are discussed explicitly, below. 


UKIRT observations were carried out with the wide-held NIR camera WFGAM, at the 
UKIRT telescope on Manna Kea, Hawaii. WFGAM uses four Rockwell Hawaii-II HgGdTe 
detectors. The held of view of each detector is 13.6' x 13.6' and the four arrays are separated 
by gaps whose areas are 94% of one detector. The observations were in H band, with each 
epoch composed of sixteen 5-second co-added dithered exposures (2 co-adds, 2 jitter points 
and 2x2 microsteps). The 18 survey helds were observed 5 times per night. 


The Wise group used the recently installed Jay Baum Rich 0.71m Telescope (G28) at 
Wise Observatory in Israel, equipped with a 1 deg^ camera. The 4 survey helds were observed 
5 times per night. At the time, the G28 had some pointing problems, and OGLE-2015-BLG- 
1285 was close to the edge of the survey held. Hence, at the initiative of ME, the Wise group 
decided to also monitor the event with the Wise Im telescope equipped with the PI camera. 
Because these observations were triggered solely to ensure coverage of an event that was in 
the survey held, with a cadence similar to that of the survey, we treat these as “survey” 
observations even though they were taken with a diherent telescope. In fact, only the PI 
observations usefully constrain the model, so we do not include the G28 data. 


KMTNet is a new survey that empl oys 4 deg^ came ras at three sites: GTIO/Ghile, 
SAAO/South Africa, and SSO/Australia flKim et al.l 120151) . While the primary goal of this 
survey is near-continuous observation of 16 deg^ in the southern bulge, it supported the 
Spitzer campaign by obtaining data on another 40deg^, with cadence of 1-2/day at each 
telescope. These lower-cadence helds included the location of OGLE-2015-BLG-1285. 


2.5. Followup Observations 

Sustained followup observations were carried out by the RoboNet team using 5 telescopes 
from the Las Gumbres Observatory Global Telescope (LGOGT) in Ghile, South Africa, and 
Australia, and by the Microlensing Follow Up Network (/iFUN) 1.3m SMARTS telescope 
at GTIO. By chance, these observations did not cover the crucial bump in the light curve, 
while their coverage of the wings adds only modestly to the survey coverage. Hence, they 
do not signihcantly inhuence the hts. They are nevertheless included for completeness. 





11 


The Spitzer team issued an anomaly alert for OGLE-2015-BLG-1285 on July 6 UT 
14:43 (JD 7210.11) based on a single OGLE data point that was posted on its web page. 
The Salerno University 0.35m telescope responded to this alert when the event rose over 
Italy just 5 hours later. These data begin at the tail end of the caustic crossing. They 
qualitatively conhrm the exit feature traced by the Wise Im. However, we do not include 
them in the £t because the target brightness was at the margin of obtaining reliable data. 

Thus, the results reported here depend overwhelmingly on survey data. 


2.6. Data Reduction 


All ground-based data were reduced using standard algorithms. Most da ta enter- 


ing the main analysis used variants of image subtraction 

Alard & Luoton 

1998 

). GTIO- 

SMARTS, UKIRT, and Wise data were reduced using DoPhot ( 

Schechter et ah 

19931. while 

the LCOGT data were processed using DanDIA Bramich ( 

20081. 


Spitzer data were reduced using a n e w algo rithm (Galchi Novati 2015, in preparation) 
whose necessity is discussed in lYee et ahl fl2015bl) . 


3. Light Curve Analysis 
3.1. Ground-based Light Curve 

The light curve contains only a single pronounced feature, which occurs near the peak 
of a roughly symmetric and otherwise low-amplitude single-lens event, i.e., less than one 
magnitude above baseline just before and after the sudden anomaly. See Figure [H In order 
to properly estimate the masses and projected separation of the components, as well as the 
distance to the system, we must unambiguously characterize the binary geometry from this 
single feature, combined with the more subtle variations of the rest of the light curve. From 
the fact that the Spitzer and ground-based lightcurves are offset by only 

Afpeak = ^peak,® ^peak,sat ~ 0.3 day, (5) 

(compared to the several week duration of the event) we already know that the microlens 
parallax effects can be ignored to hrst order for the ground-based light curve. However, 
this still leaves seven geometric parameters to be determined (to, moGe, G, «, s, g). Here, 
(to,Uo,tE) are the three parameters of the underlying single lens event, (a,s,g) are the 
three binary-lens parameters, and = pt^ is the source crossing time, where p = 6**/6*e is 
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ratio of the angular source size to the angular Einstein radius. The three underlying point 
lens parameters are, respectively, the time of closest approach to some hducial point in the 
geometry, the impact parameter (normalized to ^e), and the Einstein timescale. The three 
binary parameters are, respectively, the angle of the source trajectory relative to the binary 
axis, the projected binary separation (normalized to 6e), and the binary mass ratio. 

How can we make an exhaustive search of such a large parameter space? We begin 
by noting that the rough symmetry, dramatic outburst at peak, and lack of significant dip 
within this peak, together imply that the source is moving nearly perpendicular to the binary 
axis, and that it intercepts a cusp (or two very close and roughly parallel caustics) on this 
axis. We initially ignore the second, rather special geometry. Then, for each pair (s, q) there 
are either two cusps (for close and resonant binaries) or four cusps (for wide binaries) on 
the binary axis. Some cusps can be excluded because perpendicular trajectories would yield 
pronounced bumps as the source crossed or passed nearby to neighboring cusps. For each 
of the remaining cusps for a given {s,q), we first set the origin of the coordinate system at 
the cusp (rather than the center of mass or center of magnification, as is customary). With 
this parameterization, (to, Uq, t^) are approximately uncorrelated so that at fixed (s, q), there 
remain only two parameters with significant correlations, and a. 


We search for solutions using the Monte Carl o Markov chain technique . For each set of 
trial parameters we employ contour integration flGould fc Gaucherell 119971) for points that 
either straddle or are very close to a caustic, using 10 annuli to allow for limb darkening. For 
points that are further from th e caustics we progressively use the hexad ecapole, quadrupole, 
and monopole approximations fjPeicha &: Hevrovsky|2009l: lGouldll2008f). We u s e line ar limb- 
darkening coefficients of uj = 0.61, ur = 0.42, and M[ 3 , 6 ] = 0.28 from I Glared (120001) . based 
on the source type derived from Figure |2l In the last case, we must extrapolate. For each 
model geometry and each observatory, i, we fit for a source ffux j and a blend flux ^ that 
minimizes the of tho observed fluxes T) obs(f) relative to the predicted fluxes 


■f^i,pre(f) fs,i^(tj ^q, Uq, ^E) G; S, q^ -\- fh,i- 


( 6 ) 


We choose an initial seed trajectory by (to, mq, ^e, G, «) = (7209.75, 0, fs, 0.45 day, 90°), 
where ^e = 30 day or ^e = 50 day as discussed immediately below. For the hrst 100 trials, we 
allow only uq to vary because the crossing may be 0.02-0.1 Einstein radii from the cusp (i.e., 
Mq = 0), depending on the topology. Then, because (to, WoiG) are approximately correct and 
because there are only two correlated parameters (^e, «), the Markov chain arrives near the 
minimum very fast. 

There are only two topologies that yield a competitive x^- This agrees with the results 
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of two completely independent, and generalized search algorithms (veEI and Cl6) that do 
not make nse of the detailed featnres of the light cnrve ontlined above. In both topologies, 
the trajectory passes throngh a cnsp of the canstic strnctnre associated with the lower-mass 
component. In the hrst, it passes by the inner cnsp of the lower-mass component, i.e., the 
cusp that lies closer to the higher-mass component. In the second, it passes by either the 
outer cusp of this caustic, or essentially the same cusp of a resonant binary, i.e., the cusp 
associated with the lower-mass component. The first topology (see Figure [3] for several 
examples) typically has timescales ~ 30 days, so we use this value in our seeds for this 
topology. The second topology has ts ~ 50 days, so we seed with this value in those cases. 
(While the second topology is not shown in Figure [3l it would consist of the source passing 
roughly perpendicular to the binary axis and right through the bottom-most cusp.) 

The “inner cusp” topology is favored by = 55, which is strong evidence in its 
support. Figure m shows y^(s, q) for this topology, and Figure [3] illustrates a range of caustic 
morphologies drawn from this minimum. In addition to having different topologies, the two 
solutions are characterized by very different amounts of blend flux. As we will see in the 
next section, this implies that the inner topology is strongly favored by another, independent 
argument based on astrometric and chromatic constraints. 


3.2. Astrometric and Chromatic Constraints 

Neither the position nor the color of the apparent source changes perceptibly as it 
increases its brightness by a factor ~ 15 during the event. This is exactly what what one 
would expect for an unblended source (inner topology) but requires extreme fine tuning for 
the outer topologjfl (/;, ~ 0.68 fg). 

3.2.1. Chromatic Constraint 

Figure |2] is an (/ — H,I) CMD constructed by aligning OGLE J-band and UKIRT H- 
band data. The centroid of the clump and the “baseline object” at the location of the source 


^http://www.fisica.unisa.it/gravitationAstrophysics/RTModel/2015/RTModel.htin 
^http://astroph.chungbuk.ac.kr/ knitnet/2015.html 

'^We do not discuss this solution in detail because we consider it ruled out. However, for completeness 
we note that it has (tE, G, tte.at, tte.e, n, s, q) = (46.3 day, 0.42 day, 0.000,0.004,90.1°, 1.56, 7.2) and therefore 
{Ml, M 2 ) = (13.7,1.9) Mq. 
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are marked. We note that the if-band zero point is not fully calibrated, but such calibration 
is not needed in the present context because all results are derived from relative photometry. 

Model-independent regression of H-hand on J-band flux during the event yields 

6{I -H) = {I - H), -{I- H\ ~ (-0.009 ± 0.008) x _ ( 7 ) 

Regardless of the degree of blending, the source color is the same as that of the baseline 
object shown in Figure [21 whose position relative to the clump is 


A[(/ -H)J] = [(/ - H), J]base - [(/ - H), J]da„p = (0.02, 0.10). 
Since (/ — H)s — {I — Ff)base — 0 to high precision, this permits us to derive 


6** = 6.01 /ias 


fs 


/ba 


( 8 ) 


( 9 ) 


using the standard method of offset from the clump flYoo et ah! 120041). That is, for a 


star at the cente r of the bulge clump {{V — I),I)sfl = (1.06,14.46) (jBens 


w et 


Nataf et al.ll2013l). the sou rce size (if unblended) would be 6 ** = 6.17/las flKervella et ah 


^ 1201.4 


20041: iBessell fc Brettlll988l) . Since the baseline object is fainter by 0.1 mag, this number is 
reduced by i q-o- 2 xo i_ And since i t is redder by A(J — H) = 0.02, its surface brightness is 
lower by 4% flBoyajian et al.ll2014l) . so 2% bigger at hxed magnitude. Hence, 9^ = 6.01 /ras, 
which can then be scaled to the source flux as in Equation ([2]). 


3.2.2. Astrometric Constraint 

The source position is measured from the difference image of the event at peak magni- 
hcation relative to baseline. Since there are no stars except the source in such a difference 
image, the offset can be measured with great precision, a = 0.02 OGLE pixels (each 260 
mas). The baseline position is measured from a stack of excellent images and has a precision 
a = 0.045 pixels. Hence, the combined uncertainty in the difference between the positions of 
the source and the baseline object is a = 0.05 pixels or 13 mas. The actual difference in po¬ 
sitions is 0.05 pixels in each direction. This is consistent at the 1 a level with the hypothesis 
of an unblended source. 


3.2.3. Application to Outer Cusp Topology Solutions 

However, for the outer cusp solution {fb — 0.68/s), these constraints together imply 
that by chance another red giant of very similar color to the source and less than 1 mag 
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below the clump lies within a few tens of mas of the source. There are only four ways that 
this can happen. Either the blend is directly associated with the source (i.e., they form 
a red-giant binary), directly associated with the lens (red giant companion to the binary 
lens), a component of the lens, or the additional red giant is an unassociated star that lies 
projected within 40 mas of the source. 


The prior probability for the first option (red-giant binary source) can be evaluated 
in two steps. First, the fraction of G dwarfs (the progenitor s of bulge clump stars) wit h 
companions within 0.9 < Mcomp/^prim < 1-0 is about 3% (IDuquennov fc Mavon Il99ll) . 
On the other hand, for a 10 Gyr solar-mass solar-metallicity isochrone, the mass difference 
between 25 and 250 solar luminosities (encompassing a conservatively large range of the 
giant branch) is M(25OL0) — M(25Lq) = 1.8 x Hence, the prior probability is 

0.03 X (1.8 X 10 - 70 . 1 ) ~ 5 X 10-h 


The prior probability of the second option (red giant companion to the binary lens) is 
somewhat smaller than this, since the conditional probability for a tertiary given a few-AU 
binary is smaller than the probability of a companion given the presence of a single star. 


The prior probability of the third option depends on the details of the solution. However, 
it is a maximum if the companion mass is of order one solar mass. In this case, the above 
calculation can be applied but without the binarity factor, i.e., 1.8%. Again, this assumes 
that the solution predicts Mi ~ 1 Mq or M 2 ~ 1 Mq. Otherwise, the probability is close to 
zero. 


The probability of chance projection is even smaller. First, if the source is separated 
from the centroid by < 40mas (2cr limit), then it is separated from the blend by < 100mas 
(assuming roughly equal brightness). In this field, the density of stars that are no more than 
1 mag below the clump and within |A(J — H)\ < 0.05 of a given color (in this case, the 
source color) is lOarcmin-^. Therefore the probability of such a projection is 7r(100mas)^ x 
lOarcmin-^ ~ 10“^. In sum, the total probability is about 2% for the case that one of the 
components is about 1 M© and < 0.1% otherwise. 


3.2.4- Application to Putative Blue Lenses 

One may also apply the color constraint to the inner-cusp solutions, for which the 
blending is constrained by the £t to be small, but may not be exactly zero. In particular, 
having one or both components be main-sequence stars would in itself be consistent with 
upper limits on the blended flux. However, the color constraint implies that any such stars 
must be quite faint. For example, consider F type stars, e.g., (/ — H)q ~ 0.4. Then 
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A(/ - H) ("NJ 0.7, which implies If, — > 3.5 at 3 a conhdence. This essentially rules out 

M > 1.35 Mq main-sequence stars. 

3.3. Ground-Only Microlens Parameters 


Figure 4 shows the q) surface for the “inner topology”. The black curve shows the 
boundary between wide-binary and resonant caustic topologies flErdl fc Schneider! Il993h . 


2 ^ (1 + 

1 + q 


( 10 ) 


That is, while the minimum does he in the wide-binary caustic topology, the 1 a contour 
crosses the boundary into resonant caustics. Of course, they cannot cross very far because 
then the “neck” connecting the two formerly separate wide-binary caustics would widen, 
leading to a dip in the middle of the bump, which is not seen in Figure [TJ Hence, it has 
proved unnecessary to make an independent search of this narrow-neck resonant topology, 
mentioned at the beginning of Section [HI since it is contiguous with the inner cusp wide¬ 
binary topology. 

Table [Ushows the best-fit microlens parameters (including tte, which is discussed below) 
and their 68% confidence intervals, derived from the MCMC chain density, leading to 9-^ ~ 
0.42 mas. The surface is relatively far from parabolic, and there are non-linear correlations 
among the parameters, thus these values should not be considered as standard errors. 


3.4. Microlens Parallax Prom Spitzer Light Curve 


The microlens parallax tte of OGLE-2015-BLG-1285 is quite small and can only be 
measured because of the long baseline ~ 1.25 AU provided by Spitzer. It i s straightfor¬ 
ward to incorporate Spitzer data into the binary lens fit (e.g., IZhu et al.ll2015[) . add in two 
parameters for tte, and report the result. However, it is also important to gain a physical 
understanding of how the features of the Spitzer light curve act to constrain the parallax. 


The strongest constraining feature is the time of peak, which is Afpeak = 0.3 days earlier 
than the ground-based peak (Equation ([5])). Since the Spitzer and ground peaks are both 
due to the source crossing the binary axis, we hnd from simple geometry that 


Atpeak 

tE 


At + A/3 cot a. 


( 11 ) 


Note that this is identical to the generic Equation (|3]), but with 0 —)■ a, which follows 
from the fact that the caustic is tangent/parallel to the binary axis. Equation (ITT|l can be 
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combined with Equation ([3]) to yield 

7rE,Ar + cot a7rE,A/3 = ^ ~ 0.008, (12) 

iJ_L lE 

where the subscripts At and A/3 refer to the direction parallel and perpendicular to the 
projected separation D^, which are close to east and north, respectively. 


Then, for any given geometry that is consistent with ground-based data, the position 
of tte within the one-dimensional (1-D) space defined by Equation flT^ must be determined 
primarily by the data points on the approach to the cusp. That is, a large value of A/3 
would imply that the source crosses the binary axis substantially closer in (or farther out) 
as seen from Spitzer than Earth, leading to a somewhat different cusp-approach morphol¬ 
ogy. However, because these differences are not expected to be large, and because Spitzer 
observed with only daily cadence, we expect these constraints to be much weaker than those 
in Equation ffTTj) . Hence we expect elongated error contours in the tte plane. 


Figure 0 shows these contours. For binary lenses with ground-based data t here is of¬ 
ten a degeneracy that takes (mq) tte^w) ^ — (wq, a, tte^at) fISkowron et al.ll201ll) . which is 
sometimes dubbed the “ecliptic degeneracy” because it is exact on the ecliptic. The north 
component of tte is singled out because the ecliptic happens to run east-west for bulge fields. 
One expects this degeneracy to be preserved for Spitzer because it also lies very close to the 
ecliptic. Figured shows that this is indeed the case. 


Because M = 6*E//t7rE, and 6 *e varies very little between viable solutions, the total mass 
M tracks (inversely) the tte values in Figure [5] very well. However, because a range of mass 
ratios q are permitted (see Figure H]), the probability contours for the two component masses. 
Ml and M 2 , (Figure [H]) are less 1-D than the tte contours. These contours were calculated 



Finally, we note that from the dehnition of 0 e (Equation ([T])), and the fact that ^e is 
essentially the same in all solutions, we have 

f 0E M V' 

Since the source is nearly at the center of the clump, it is almost certainly in the bulge, with 
distance Dg ~ 8kpc. Equation flT^ then implies Dg — ~ 0.6 kpc(M/2 Mq)“^, i.e., that 

the lens is also in the bulge. 
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Taking account of the ~ 8% error in 6^, (and so ^e), there is an additional 1 a error in 
M and TTrei of 8% that must be added in quadrature. However, this is smaller than the errors 
that propagate directly from htting the light curve. 

Figure [7] shows the probability distribution of the primary mass Mi, which is peaked 
at Ml ~ 2.0 Mq (black). The fraction of the curve area M > 1.35 Mq (typical NS mass) is 
80%. As mentioned in Section 13.2.41 main-sequence stars in this mass range are ruled out 
by the chromatic constraints. Together, these imply that this primary is a massive-remnant 
(NS or BH) candidate. 

Table [2] summarises the median and the 68% conhdence intervals for the physical param¬ 
eters of the binary system. As for Tabled! we warn about the non-gaussianity and non-linear 
correlations among the MCMC variables, which are also reflected in the uncertainties on the 
physical parameters. 

4. Future Mass Determination from Proper Motion Measurement 


The nature of OGLE-2015-BLG-1285La, i.e., whether it is a massive remnant and if so 
whether it is a BH or NS, can be decisively resolved by a proper-motion measurement of 
its companion, OGLE-2015-BLG-1285Lb, whose mass implies that it is almost certainly a 
luminous low-mass main-sequence star or a white dwarf (see EigureEj). This is illustrated 
in Figure O where we show a hypothetical future measurement of the direction of lens- 
source relative proper motion /x with an error of either 3° of 1°. Because vTrei is very small, 
/X = AO/At, where AO is the observed lens-source separation at a future epoch tpeak + Ah 


In order to make such a measurement, it is of course necessary for the source and lens 
to separate. Since their proper motion is known u = Ov./t-p, ^ 4 . 8mas yr~^, the wait time 
depends primarily on the resolution of the telescope. iBatista et al.l (120151) were able to resolve 
the source and lens of OGLE-2005-BLG-169 using Keck when they were separated by only 
61 mas. However, in that case the source and lens had comparable brightness, whereas here 
they have a flux ratio fi/fs ^ 1%- Therefore, it is likely that 15 years would be required with 
present instruments. However, during this interval, it is likely that 30m class telescopes with 
AO capability will come on line. By the time that they do, this measurement will already 
be quite feasible. For example, for separations of A9 ~ 50 mas and limiting resolution 
FWHM ~ 11 mas (e.g., in J-band on the Giant Magellan Telescope), a 1° measurement 
would require that the companion be centroided to 8% of the FWHM in the transverse 
direction. Improvements to higher precision would be considerably more difficult because 
the orbital motion of the secondary about the center of mass is of order 0.5 mas. 


Figure [7| shows the impact of such a proper motion measurement on the estimate of the 
primary mass for ±3° (brown) and ±1° (magenta) errors, respectively. Such a measurement 
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would automatically detect (or rule out) light from a main-sequence (or giant-branch) pri¬ 
mary. If the primary is indeed dark, then the detection of light from the secondary would 
yield a mass estimate of that body, which would in principle constrain the solution and so 
further constrain the mass of the primary. However, in practice we hnd that this does not 
lead to signihcant further improvement beyond what can be achieved with the proper motion 
measurement alone. 


Finally, we note that iGouldl (1201411 has investigated the problem of turning 1-D parallax 
measurements into 2-D parallaxes via such proper motion measurements and shown that, 
in general, there is a discrete degeneracy induced by the difference between the geocentric 
frame of the tte measurement and the heliocentric frame of the measurement. While 
the physical origins of the 1-D degeneracy are completely different in the present case, the 
mathematics, leading to a quadratic equation, are identical. However, in t he language of 
that paper, H <C 1 and (7 -C 1, so the “alternate” solution (Equation (10) of iGouldl 120141) is 
so large as to be easily ruled out by the light curve. Hence, there is no degeneracy. 


5. The General Problem of Breaking the 1-D Binary-Lens Parallax 

Degeneracy 


As discussed in Section [H events that have only a single caustic crossing that is mon¬ 
itored from both Earth and a satellite generically suffer from a 1-D degeneracy in tte. See 
Equations dH) and flTT|) . Breaking this degeneracy, in particular breaking it at the high preci¬ 
sion required for bulge-BH mass measurements (with their very small tte), requires additional 
high-precision information. In Section dl]) we showed that for the specihc case of OGLE- 
2015-BLG-1285, it will be possible to break this degeneracy with future proper motion mea¬ 
surements, although this still will not deliver a precise mass measurement. Moreover, Erst, 
this requires waiting many years, and second it will only be possible because the secondary 
(or possibly primary) is luminous. Hence, this method will not be applicable to binaries 
composed of two remnants, which is the unique province of microlensing (and arguably the 
most interesting case). 


Thus, it is of interest to understand how this degeneracy can be broken from microlensing 
data alone, and in particular, what can be done to modify current experimental protocols 
to increase the chances of success. Generically, this requires an additional feature in the 
light curve that is monitored from both Earth and the satellite and whose appearance from 
the satellite is predicted to be a function of the Earth-satellite separation vector projected 
onto the Einstein ring. For example, wh ile OGLE-2014-BLG-1050 does suffer from a 1-D 
degeneracy (Figure 3a of IZhu et al.l 1201511 . it is much less severe than that of OGLE-2015- 
BLG-1285. This is primarily becaus e the source trajectory skirts the interior edge of the 
caustic (Figure 2a of IZhu et al.ll2015l) . which has signihcant structure. Another widespread 
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feature of binary light curves that could provide such leverage is a cusp-approach “bump”, 
which often occurs shortly after a caustic exit. 


However, by far the best structure would be a second caustic crossing. In general, if 
there is a caustic entrance, there must be a caustic exit. This was not actually true of 
OGLE-2015-BLG-1285 because the source was much bigger than the separation between the 
caustics, so the entrance and exit combined to form a single bump. However, if the source had 
been a main-sequence star, i.e., 10 times smaller, there would have been a dip between the 
the entrance and exit, whose duration would probe tte in the orthogonal direction (although 
not with our data set, which has only two points over the entire cusp crossing). 


More generally however, if there are two independent crossings with measured Afi 2 and 
well-determined crossing angles 0 i ^2 (derived from the model), then one easily hnds 


I cot 01 — cot 02 1 I cot 01 — cot 02 1 

where we have assumed that the errors a{Ati^ 2 ) in the determinations of Afi ^2 are the same. 
The most challenging case is clearly |A0| = |0i — 02| 1, for which 


cr(Ar) 


sin(20) a{At) 
\/2|A0| ^E 


cr(A0) ^ 


1 — cos(20) a{At) 
\/2|A0| 


(15) 


It is quite plausible to reach errors of ~ 0.01 days for these time offsets, so that even assuming 
that the hrst factors in Equations flTdl) and ffTbll are of order ~ 10, the errors cr(Ar) and 
ct(A 0) would be only ~ 0.002 for a ^e ~ 50 day event. Thus, coverage of two independent 
caustic crossings is by far the best method to measure tte for binaries. We discuss how this 
can be achieved in practice in the next section. 


6. Discussion 


While it is not yet known whether OGLE-2015-BLG-1285La,b contains a massive rem¬ 
nant, its example at least shows that such detections are possible. The main challenge to 
detecting more of these systems is simply to monitor a large number of targets, including 
both light curves that already have clear binary signatures, but also those (like OGLE- 
2015-BLG-1285) that erupt with these signatures unexpectedly. Spitzer has the advantage 
that any known microlensing event from the ~ 100 deg^ that are currently monitored can 
be targeted. In contrast to other classes of interesting microlensing events, most particu¬ 
larly free-floating planets, events containing massive remnants tend to be quite long, so that 
neither the low survey cadence in the majority of these helds nor the relatively long time 
required to upload targets presents a serious obstacle. A favorable feature of these binaries 
relative to other microlensing binaries is that, due to their typically small micro lens parallax 
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tte = \/vTi-ei/ kM, if the binary is seen to be in a caustic trough from Earth, it is likely to 
also be in the trough as seen from Spitzer. Hence, it is likely that at least one of the sharp 
features induced by caustic crossings will be monitored. However, as discussed in Section |5l 
such events will generally not lead to a precise mass measurement unless there are additional 
features in the light curve, such as a post-caustic-exit bump due to a cusp approach. This 
implies that systems that are monitored from before the caustic entrance, and therefore not 
usually known to be binaries (e.g., OGLE-2015-BLG-1285) will be the most favorable for 
precise mass measurements. 

Such measurements for these systems would greatly beneht from more aggressive ob¬ 
servations from both the ground and Spitzer. In the present case, we were fortunate that 
the microlensed source was a clump giant, with roughly 10 times larger radius than the Sun. 
This meant that the caustic features lasted 10 times longer than for a solar-type source. 
Given the relatively sparse coverage over the peak, it is possible that critical portions of 
these features would have been missed entirely from the ground if the source had been sim¬ 
ilar to the Sun. The Spitzer observations were even sparser, with only two over the bump. 
There are a number of modihcations that could be made to more than double the cadence 
of Spitzer observations (given the same overall time allocation). In particular, with new 
real-time reductions (Galchi Novati 2015, in prep), it should be possible to stop observations 
of many events whose Spitzer light curve has essentially reached baseline. In addition, it is 
probably more productive to increase overall cadences at the expense of the 2015 campaign’s 
extra observations for events that are relatively high magnihcation as seen from the ground. 

For OGLE-2015-BLG-1285, surveys were almost entirely responsible for capturing the 
peak of the event from the ground (see Section [2]). In this respect, two surveys that were 
specihcally created for the Spitzer campaign (UKIRT and Wise northern bulge surveys) as 
well as the wide-area low-cadence survey undertaken by KMTNet in support of this program, 
played an important role. The continuation of such surveys and the organization of new ones 
will be crucial. In particular, we note that the VISTA telescope is well placed to do such a 
survey in highly extincted bulge regions. 

However, it is also the case that followup observations could greatly enhance the chances 
of detecting caustic features in areas that are not well covered by surveys. The primary 
motivation for such followup observations to date has been planet detection. This both 
drives the allocation of available followup resources, which are focused on planet sensitivity 
and detection, and also fundamentally limits the total amount of followup resources to 
those available to planet hunters. For example, the very intensive LGOGT Spitzer support 
campaign ran out of allocated observing time (due to “too much” good weather!) two days 
before the peak of OGLE-2015-BLG-1285. While the LGOGT team did arrange to get some 
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additional points of this event in response to the anomaly alert, the main point here is that 
detection and characterization of BHs is very challenging and observational resonrces have 
been limited partly becanse the potential to detect BHs is not widely appreciated. 

The K2 microlensing campaign, schednled for 83 days beginning April 2016 will pro¬ 
vide a nnique opportnnity for space-based microlensing withont the need for gronnd-based 
alerts, and hence with a greater chance that canstic entrances will be monitored from space. 
This advantage (relative to Spitzer) is balanced by the fact that these events will be drawn 
from a relatively small area, albeit one with close to a peak snrface density of microlens¬ 
ing events. From the standpoint of making binary-lens mass measnrements, and BH-binary 
mass measnrements in particular, we note that it is exceptionally important that this entire 
area be monitored from the ground at high cadence and as continuously as possible. For 
example, even extremely faint stars can give rise to briefly bright caustic crossings that can 
be effectively monitored by Kepler with its 30 minute cadence, even if the majority of the 
light curve cannot. However, only if the corresponding caustic crossings are monitored from 
the ground well enough to effectively model the light curve, will this result in accurate mass 
measurements. A very aggressive attitude toward continuous coverage will be especially im¬ 
portant toward the beginning of the campaign when individual southern sites can observe 
the bulge for only hve hours per night. 

Gaia will provide complementary information on BH binaries. There are ~ 3 x 10® 
G dwarfs within 250 pc of the Sun, and it is hoped that these will all have a{7i) ~ 7/ras 
parallaxes by the end of the 5-yr mission. For those that have BH companions with periods 
less than the mission lifetime, the semi-major axes of their orbit about the binary center of 
mass could be measured with the same precision. This would imply 7 a detections for all 
those in the semi-major axis range 0.015 AU < a < 5.3 AU (assuming Mbh = 5Mq). Even 
the 1.5-year data release would enable detection of those in the range 0.027 AU < a < 2.4 AU. 


Of course, Gaia cannot detect systems that, like OGLE-2015-BLG-1285L, lie in the 
Galactic bulge, nor can it detect totally dark systems to which microlensing is sensitive. 
However, WFIRST will enable two different probes of BHs that explore both of these regimes. 
First, it will be able to hnd and measure the mass of isolated BHs, as well as BH binaries by 
combining microlens parallax measurement s from its superb p hotometry with astrometric 
microlensing from its excellent astrometry ( Gould fc Yed 20141). For discu s sion o n possible 


astrometric microlensing measurments with Gaia see iBelokurov fc EvansI (1200211 . Second, 
i t will obtain ctItt) < 4 pas astrometry on of order 4 x 10^ stars in its microlensing helds 
(iGould et al.l 1201511 . The majority of these stars will be in the Galactic bulge, where this 
precision corresponds to a 7 a threshold of ~ 0.25 AU. Hence, WFIRST will be sensitive to 
BH companions to ~ 4 x 10^ luminous stars over the range 0.3 AU < a < 5.3 AU. 
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Table 1: Best-fit microlensing model parameters and their 68% uncertainty range derived 
from the MCMC chain density (for both mq > 0 and Uq < 0). We note there are non¬ 
linear correlations among the parameters, and that the entire surface is not parabolic. 
The solutions are almost symmetric, with a small difference in tte, due to the small offset 
between the projected Spitzer-Eaith axis and the ecliptic. 


Parameter 

Mo > 0 

Mo < 0 

to - 2457200 [d] 

9.74 [9.73,9.75] 

9.74 [9.73,9.75] 

Mo 

0.46 [0.42,0.53] 

-0.46 [-0.42,-0.53] 

[d] 

31.4 [29.5,32.5] 

31.4 [29.5,32.5] 

TTE.N 

0.019 [0.012,0.029] 

-0.018 [-0.012,-0.029] 

tJ'e.e 

0.0087 [0.0075,0.0112] 

0.0092 [0.0080,0.0121] 

t* [d] 

0.455 [0.450,0.459] 

0.455 [0.450,0.459] 

<1 

2.9 [2.1,3.4] 

2.9 [2.1,3.4] 

s 

1.934 [1.928,1.944] 

1.934 [1.928,1.944] 

a [deg] 

80 [78,82] 

-80 [-78,-82] 


Table 2: Physical properties of the binary system. Median and 68% uncertainty range values 
derived from the MCMC chain density, after accounting for the Jacobian of the transforma¬ 
tion between the MCMC variables and the physical quantities. 


Parameter 

Median 

68% conhdence intervals 

Ml [Mq] 

2.0 

[1.2,3.3] 

M2 [Mq] 

0.8 

[0.5,1.2] 

[AU] 

6.1 

[5.7,6.5] 

Dl [kpc] 

7.5 

[7.3,7.7] 
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Fig. 1.— Light curve of OGLE-2015-BLG-1285 with data from Spitzer (red) and various 
ground-based observatories (see interior hgure labels). The magnitude scale is what is di¬ 
rectly observed by OGLE. All other observatories (including Spitzer) are aligned so that 
equal “magnitude” reflects equal magnihcation. The very small (~ 0.3 day) offset between 
the peak as seen by Spitzer and the ground hints that the microlens parallax tte = 
may be small, which would imply a high-mass lens. For ground-based data, two models are 
shown, one for if-band limb darkening (blue), which should be compared to UKIRT data, 
and one for /-band (black), which should be compared to all other data. The difference 
between those two curves can be seen only at the peak of the anomaly. 
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(l-H) [instrumental] 


Fig. 2.— CMD derived by combining calibrated /-band photometry from OGLE with 
instrumental //-band photometry from UKIRT. The centroid of the clump (red) and 
the “baseline object” (blue) are marked. From this offset one derives the source radius 
6** = 6.01 p-as(/s//base)^^^, where fs/fhase is the ratio of source flux to baseline flux. 
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Fig. 3.— Three different caustic structures that are consistent with the light curve (Fig. [T]), 
i.e., from the 2 a region of Figure HI All three are rotated by a and scaled by so that the 
x-axis is simply source position as a function of time. Hence, the binary axis is oriented so 
that the primary is toward the top of the plot and the secondary is toward the bottom. Both 
wide-binary and resonant caustic topologies are permitted. Main panel shows full caustics 
with Spffeer trajectory shown for one of the three caustic structures. Inset is a zoom showing 
the source size and its position at times of observations, with same color scheme as Figure [H 
in which Spitzer points are shown relative to the caustic rather than clock time. An alternate 
topology, which is ruled out by arguments given in Sections [3.1 l and lT^ would have the source 
trajectory pass through the bottom-most cusp, nearly perpendicular to the binary axis (not 
shown). 
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Fig. 4.— surface for “inner cusp” topology, where s is the projected separation of 

the components normalized to 6e and q is the mass ratio of the primary to its companion. 
Black curve = (1 + + q) shows boundary between wide-binary topology (two 

4-sided caustics) and resonant topology (one 6-sided caustic). See Figure O 
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7rE(East) 

Fig. 5.— X^('^e) surface. Solutions with mq > 0 (mq < 0) lie in the upper (lower) part of the 
diagram. Physical explanation for quasi-l-D contours is given by Equations flTT]) and flT^ 
in Section 13.41 The magenta and brown rays converging at the origin show the 
impacts for hypothetical future proper-motion measurements fi = (60 ± 1 )° and 
jl = (60 ± 3 )°, respectively. 
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Mi[Mq] 


Fig. 6.— Binary mass map. Contours showing 68%, 95% and 99.7% probability regions. 
The most probable combination is a NS primary and a main-sequence secondary. However, 
binaries with two main-sequence stars or two massive remnants (e.g., BH-I-NS) are also 
possible. 
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Ml [Mo] 


Fig. 7.— Differential probability distribution of the primary mass Mi, assuming a prior 
uniform in log(M), as of “now” (black), i.e., based solely on microlensing measurement. 
About 80% of the probability lies above Mi > 1.35 Mq, making this a massive-remnant (NS 
or BH) candidate. Magenta and brown curves show the impact of future proper motion 
measurements jl = (60 ± 1)° and jl = (60 ± 3)°, respectively, as indicated in Figure [5l 










